TRANSPORT PARAMETERS 
FOR THERMAL 
NEUTRONS IN WATER 

BY 

FRANCIS JOSEPH SISK 



S 54* 



Library 

V. 8. Naval Postgraduate School 
Annqpoks, MaL 



TRANSPORT PARAMETERS FOR THERMAL 
NEUTRONS IN WATER 



by 

Francis Joseph Sisk 
Lieutenant, United States Navy- 



Submitted in partial fulfillment 
of the requirements 
for the degree of 
MASTER OF SCIENCE 



United States Naval Postgraduate School 
Annapolis, Maryland 
1950 



:KRTiriCATE OF APPROVAL 



This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 



from the 

United States Naval Postgraduate School 



TABLE OF CONTENTS 



Page 



Chapter I. Introduction and Summary 1 

1. Introduction 1 

2 . Summary 4 

Chapter II. Transport Theory and Methods of Analysis 5 

Chapter III. Experimental Procedure 22 

Bibliography 32 

Appendix I 33 

Appendix II 37 



(iii) 



LIST OF ILLUSTRATIONS 



Page 

1. Velocity Diagram for Two-Body Collision 26 

2. Integration for Plane Kernel - Transport Equations 27 

3. Foil Mount and Cadmium Plane Assembly 28 

4. Thermal Neutron Flux in Water Near a Black Plane 29 

5. One -Dimensional Thermal Neutron Flux in Water Near 30 

a Black Boundary 

6. Diametral Flux Distribution 31 



(iv) 



CHAPTER I 



INTRODUCTION AND SUMMARY 

1. Introduction 

The steady state spatial distribution of neutron fluxes in matter is a 
problem of prime importance to the designer of chain reacting piles, and 
has been the subject of very much theoretical and experimental work. The 
case in which the neutron energy is large compared to that of the nuclei 
constituting the medium (the "slowing down" case) is of first importance 
and is rather extensively investigated. The alternate case in which the 
neutrons are in energy equilibrium with the medium is of lesser importance 
and comparatively little experimental work has been done, although the 
fundamental theory is well in hand. The cross-sections for neutron scat- 
tering and absorption appropriate to the first case are the "ordinary" 
energy-dependent cross-sections (or simple functions of these) which have 
been so extensively investigated by various neutron spectrometry groups* 

The cross-sections appropriate to the second case are more vaguely defined. 
Perhaps as adequate a definition of these parameters as any is the "ad hoc" 
one - those numbers which, when inserted in the mathematical theory of 
particle transport yield proper solutions for the neutron flux. It can 
be demonstrated that these numbers are not in general the average values 
of the energy dependent cross sections weighted by the appropriate energy 
distribution function for the ambient temperature. 

The "ad hoc" definition suggests a means of experimentally measuring 
the numbers. One requires a medium - of simple geometry so that the trans- 
port problem is soluble - in which a flux of thermal equilibrium neutrons 
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is maintained and in which the flux can be measured at various points with 
neutron detectors o Since two parameters are to be determined, two independent 
solutions for the flux must be available - but this is easily provided for 
by considering (l) a region "remote" from discontinuities in which the 
solution is asymptotic and (2) a region adjacent to a discontinuity - say 
a "black" boundary . One measures the flux at a sufficient number of points 
in each region to define the spatial functions. The transport problem is 
solved many times for various values of the parameters, and if fortune is 
good one pair of values will be found to yield solutions which fit the ex- 
perimental functions in both regions. 

To particularize = the chosen medium is normal water, the geometry 
one -dimensional consisting of a tank many mean free paths in its dimensions, 
traversed by a uniform thermal neutron current arriving at one face and 
absorbed by a layer of cadmium at the opposite face. The neutron flux is 
measured at points along a line normal to the cadmium facing and extending 
from the facing well into region (l). 

The one-dimensional transport problem is integrable in closed form for 
region (l). This solution is matched to the experimental data, yielding a 
characteristic equation in say the total cross-section and the scattering 
cross-section - thus the scattering section may be eliminated. The trans- 
port problem is then integrated numerically in region (2) for various 
values of total cross section and that value selected which yields best fit 
with experimental data. Unfortunately, this straight-forward procedure 
could not be carried out. It was found that numerical integration produced 
a function which fit the experimental data astonishingly well, but that the 
values of the parameter which yielded this good fit were variable with 
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respect to the integrating procedure. Further, it was found that each 
improvement in the integrating plan led monotonously to a reduction in 
the total cross section, indicating that the process of improving the 
operation was not near convergence. The rapid increase in computing labor 
made it necessary to abandon this process before a firm solution was ob- 
tained. 

The Wiener-Hopf method, treated by Frankel and Goldberg (3), offers 
a general solution to the one-dimensional transport problem in a single 
semi-infinite medium. However, it is available only as an exceedingly 
difficult double numerical quadrature which could not be handled for lack 
of computing facility. A solution is given for the flux at the black 
boundary only, requiring only one numerical quadrature. By equating this 
solution in terms of the parameters to the experimental value of the 
"contact" flux, unique values are obtained for the parameters. These 
values extrapolate the trend of the numerical solutions, and are in that 
sense to be regarded as supported by the numerical method. 

The convergence on the Wiener-Hopf solution may be illustrated by 
the successive values for the macroscopic "total" cross section. They 
ares Numerical - 2.9 cm"*) 2.5 cm”* and 2.3 cm"*j Wiener-Hopf - 2.10 cm“* . 
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2. Summary 



Findings may be tabulated as follows: 



Parameter - Water at 24° C. 


Values 


Diffusion Length 


2.67 ± .02 cm. 


Transport Mean Free Path 


0.48 ± .01 cm. 


Absorption Mean Path (Equivalent) 


44.2 * 1.3 cm. 


Microscopic Transport Cross Section 


62.4 ±1.5 barns/mol. 


Microscopic Absorption Cross 
Section (Equivalent) 


0.68 0.01 barns/mol. 


Ratio - Boundary Flux -r Asymptotic 
(Exponential) Flux 


0.186 


Ratio - Boundary Flux 4- Asymptotic 
(Hyperbolic) Flux 


0.87 


Extrapolation Distance = X 0 


0.32 cm. 


Ratio - Extrapolation Distance/Total 
Mean Path 


0.67 



Comparison of the cross-section with the energy-dependent values for 
protons in water by Bacher(l) and Melkonian (7) indicate an "average" 
neutron energy for transport of perhaps 0.05 electron volts in contrast 
with the Maxwell-Bolt zmann mode of 0.025 e.v. This is not surprising since 
it is expected that transport is largely effected by the high energy neutrons 
which have longer free paths. In addition it must be remembered that the 
neutron flux is weighted by a velocity term so that the average flux corres- 
ponds to a higher energy than does the average neutron, and transport weights 
the cross-section by the flux rather than by the neutron density. The l/v 
absorption in hydrogen provides a preferential sink for low energy neutrons, 
thus the thermal distribution is displaced from the Maxwell-Bolt zmann toward 
higher energy. 
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CHAPTER II 



TRANSPORT THEORY AND METHODS OF ANALYSIS 

To indicate the assumptions involved and to portray the physical model 
which was analyzed, it may be well to review the fundamental theory of 
neutron transport. It must be emphasized that the model of neutron trans- 
port for neutrons in energy equilibrium with the medium, at energies well 
below the atomic binding, differs radically from that for energetic neutrons 
(E n > 1 ev) in media at normal temperature. 

In the latter case the energy transfer in neutron -nucleus interactions 
is uniformly from neutron to nucleus, the nuclei interact as free particles, 
and since the atomic vibrational energy is comparatively small, the nuclei 
may be regarded as initially at rest. None of these statements hold for 
the former case. 

In the present problem, by definition, the net average energy transfer 
between neutrons and medium is very small and due solely to the preferential 
absorption of the slower neutrons. The protons are bound in the molecule 
with energy larger than the neutron energy and the nuclei cannot interact 
as free particles. The thermal motion of the atoms and of the systems in 
which they are bound is not negligible, and so affects not only the mechanics 
of the collisions but also the frequency of occurrence since the collision 
probability is a function of the flux of nuclei as well as of neutrons. 

Inelastic interaction of neutrons with rotational and vibrational 
levels of the water molecule are probably present. The rotational levels 
are quite numerous and begin at wave number 23.8 cm”^ , (about .003 electron 
volts) according to Wu (9). 

Quantitive discussion of the scattering mechanism in a system as 
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complex as this is not in the scope of this paper. However to avoid mis- 
interpretation of the results, it is necessary, to consider at least the 
question of whether the scattering is symmetric, that is whether all values 



of the angle of deflection of the neutron in the Laboratory System are 
equally probable, when the neutrons are in energy equilibrium with the 
medium. 

Referring to Fig. (l), and v are the velocities of neutron and 

molecule respectively before collision, in the laboratory system, and a> 
the angle between them. To transfer to the center of mass system a vector 
velocity V is added of such magnitude and direction as to reduce the total 
momentum to zero. Using unity as neutron mass and M as molecule mass we 
have : 

v D -I- V x + M(v cos w + V x ) =0 



and ^ 


4- M(v 


sin cj 4- Vjj) 


= 0 




thus w _ 


M 


1 






v* = - 


M + l v 


cos w ~ 


r v ° 




and \, _ 


_M_ v 


sin cj 






V 3 ' 


M+l v 






where and y are the scalar components 


of V • 




The transformation 


vector is 


thus seen to resolve into two vectors, 


one of 


magnitude — — — v 

M + 1 


and 


opposite to v 


in direction, and one 


of rag- 



nitude I v and opposite to in direction. The vectors vj,' 

- M + l 

and v' are the velocities before collision in 0.; . syste/ and and V 

are the C.M. velocities after collision. If the scatterin spherically 
symmetric in C. id. system, all values of cosy, the deflection angle, are 
eq , a‘’' , ' J L ’obahlr. To "^air ’.he neu + ror velocity (J , collision, in 
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Laboratory system, we add (— V) vectorially to %“ . The probability of 
oj is sin u> , and is therefore symmetric, so on the average the contribu- 
tion of — — — v involves no preferred direction. However, the contribu- 

M + l 

tion of - — ! y inserts the direction of the incoming neutron as a 

M + I ° 

component of u . Thus the scattered neutrons will emerge from collisions 
preferentially in the direction of their pre-collision motion with a weight 
increasing with neutron energy. 

The preferentially forward scattering, as indicated, correlates rather 
poorly the direction of each path with the preceding one, less well with 
the next but one preceding, and after a large number of intervening scat- 
terings the correlation vanishes . One may attempt to obtain the number of 
neutrons passing through a small volume increment (dY) by summing up those 
neutrons scattered in all other volume increments ( d y') in the space which 
reach the first increment in free flight. However, since the number of 
scattered neutrons leaving dY in the direction of d Y depends on the aver- 
age direction of the neutrons entering dY', this average direction must be 
found by another summing up of previously scattered neutrons and so on. 

The mathematical analysis of such a situation yields infinite order multiple 
integrals, so simplification is necessary. The difficulty does not occur if 
the scattering is isotropic, so an obvious suggestion might be to approxi- 
mate the system with a model in which the scattering is isotropic. 

If one considers a neutron scattered isotropically at a point in a non- 
absorbing medium and asks for the probability of finding it at any point in 
the space after the lapse of a long time, it is clear that the resulting 
probability field is symmetric about the initial point. Analogously, one 
considers a neutron with initial motion % scattered at a point anisotrop- 
ically in the Lab system but isotropically in C.M. system. 
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The probability of finding the once scattered neutron at position r 
is then symmetric about the moving center of mass whose velocity is 

V 0 . However, the probability of finding the second scattered 

M + l 

neutron at 7 is only azimuthally symmetric since each spherical shell of 
second scattered neutron probability is immediately again retarded by the 
factor ! , corresponding to the motion of the center of mass of the 

M + l 

second scattering system. Thus the second scattered neutrons which have 
made long flight after first scattering will be found on the average farther 
along in the direction of V 0 than those which make short flights. In 
spite of the forbidding complexity of the picture some information can be 
obtained about the probability distribution after a large number of collisions. 
The path lengths X , the azimuthal scattering angles JL , and the cosines 
of the polar scattering angles jx , are every one statistically independent, 
and may be averaged independently. The centroid of the distribution is 
obviously found in the direction of v 0 and may be located by first averaging 
over the azimuthal angles starting with the last flight. One obtains 

+ Ki fifx + ••• >?n II, /**+•••) * Averaging over i n and 
one has 



7 ^ d /*» ^2 d *2 + * * 



The required first moment is then > J 

which is the expansion of . One concludes then that the average 

l_ ^ 

position, after very many scatterings, of a neutron projected from a point 

in a medium is at a distance A s (l + , ) or 7 ^ 5 _ from the 

I I “/A 

initial point. This parameter is called the transport mean free path. The 
higher moments of the probability function are not available in general. 
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F. Joliot (5) has computed the probability function itself for the nth 
isotropic scattering in a quite complicated series of Hermite polynomials. 

One may now reasonably postulate that the scattering mechanism in 
an anistropically scattering medium may be rather well reproduced by an 
isotropic model, provided that the scattering length be weighted by the 
factor 1 o This ignores the defect that the higher moments are not 

1 -F 

properly scaled. 

It must be emphasized that capture has not been considered and that 
the model is not correct for large capture. It then remains to be demon- 
strated that the model is satisfactory if capture is small but not negligible. 

The simplified transport model now may be more clearly defined. It 
consists of a three dimensional array of fixed lattice points at which the 
neutrons scatter isotropically, or are absorbed. 

The scattering length in the model is the transport free path in the 
real medium and is in general not known. The absorption free path in the 
model is not in general the same as the absorption length in the medium but 
is probably nearly the same if absorption is small. A more useful pair 
of parameters are the inverses of the free paths - the macroscopic cross 
sections - defined as the probability per centimeter of path of the defined 
event. To prove the relations 



“ d 



I = 
I = 



I 2 cU 



Iofc 



-2 J? 



Let I 0 = I } 

The mean value of 



I - t' 

J( is s 



-ijt 



Where I is the beam intensity 



r*-« d; 



I 



9 



A macroscopic absorption cross-sectior is similarly defined. 

The following notation is to be used: 
n - Neutron density, cm"® 

- Neutron scalar flux density, called simply "flux,” cm"^ sec - ^ 

- Mean macroscopic transport cross section, cm"^ - 

- Mean macroscopic absorption cross section cm"'*' 

- Mean total cross section, cm - ^ 

- Lattice point density of the medium, cm"® 

<T - Microscopic total cross section, cm^ 

% - Distance, cm. 

£ - Distance in mean paths 

F - Position vector 



It 

la 

L 

N 



One asks for the number of neutrons scattered in incremental volume 
at position F and crossing unit normal area at r 7 per second. The total 
scattered at r is ft(r)E±dY - all cross sections being independent 
of F in the homogeneous medium - and of these — dco leave r in d w 

4 n 

about solid angle u> . Since the area subtended at F 7 by dw at r is 
(F 7 — r ) 2 doJ , a division by this factor is indicated, and since the 
beam in any dw is attenuated by absorption and re-scattering as 



9 multiplication by this last factor 

is required. The required number is then -Z P-r 

V?(r) 2 . t fe 1 1 01 r 

4* (r'-r) z 

The neutron flux at r 7 is then obtained by summing up the neutrons 
scattered at all r and crossing unit area at r . Thus 



^ f 

Am 






-l|r'-r| 



(r'“ r ) l 



dY 
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In the problem at hand, the medium is considered to be a semi-infinite 



slab - infinite in horizontal extent and downward, its upper face bounded 
by a plane from which no neutrons return, and having a finite neutron current 
passing upward through it. The coordinate % is measured from the "black" 
plane and normal to it. In this system the flux is a function of X alone. 
One requires the number of neutrons scattered in an incremental slab 
dx at x and passing in free flight through unit area about a point at % / . 
Referring to Fig. 2, this quantity called is found as follows* 




r 2 + (%'- xf = p 2 



r dr = p dp 



K 




K 
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The integral is termed E, » the exponential integral, 

and is identical with the function as tabulated by Jahnke & Ehmde 

and others . 

The flux at x' is obtained by summing the contributions of all increments 
dx at % as follows: 

= -^-J <P(x.) E,( S | X '-z|)dx 

"space 

This linear integral equation is perhaps more conveniently handled in mean 
free path dimensions as follows: 



%) = f (P(i) E i (ii-.ii) d 2 



For the configuration space described, the limits of integration are zero 
and plus infinity. However as -Z — «-■'* 15 the contribution to the 

integral between negative infinity and zero becomes negligibly small, and 
the limits may be taken, with vanishing error, over all values of Z 
In this asymptotic form the equation has a simple hyperbolic solution. 

The usual method of treatment, to indicate the form of the asymptotic 
solution (as in Frankel & Nelson (4)), involves transforming the origin to 
write : 

Vp) m rtf <V+u 

-'-oo 



and expand in Taylor's series about H S 

^'1 = ft f E ' 0*1) + ZPp'l + f[ ?"(!') + -]di 

<30 



TD 



/ oo 

E 



i (i*i) 



da 



is the moment of the kernel 
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about its origin, and this vanishes by symmetry for m odd. One has then 

*V) = *Vd 

neglecting terms above the third order. Rearranging, one has the one- 
dimensional differential equation of elementary diffusion theory - the 
Helmholtz equations 

2TQ 0 _ 41 



m " _u ( \ u) = o 

" + \ rrix 



for which the solutions are hyperbolic if IT) 0 ^ 
7f| o > 2^— , and linear if TT) 0 = ^ 



21 



circular if 
Now 



m c 



r°° 

= J E l(|?|) d2 = 2 , and if absorption is present, then , 



. Thus we expect hyperbolic solutions, so i 



s ° 7n 0 = 2 < 

t _ 2 

trial function ^ is substituted in the integral equation, and 

the integration performed at the point Z /= ° 



I = 



It 

21 



It 

21 



It 

21 



-tfz , 



dz 



'* E| (ls'-il) 

'{//** E, f-« da 



-{f 






'(*) 



dz- + 






~^Z~f c ° sh (* ? ) E .(z) dz 

'0 

by parts after putting (?) in integral form. 



which is integrated 



, 



/ oo /-CO _ u 

cosh(JCz)ds j — - cl ^ 

J?- y 
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is of 



integral vanishes for J^< | since at upper limit £»(f) 



and at lower limit 



sinh 

a 



goes as 2- and the singularity of 



£ | is logarithmic. Writing the second integral in exponential form 



= _£t_ 
2*Z 



f 

-'O 



t -g(i-x)_ t -z(n-*) d! 



and separating the range of integration into intervals: 



I • kr\l 4 

1° 0 »-<* *nfer J 

the third and fourth integrals are seen to he each E]^(l) and cancel, thus 



Zt o l + <# 
2 /£ l-tf 



It - 

Z 






expanding in power series 

z 



\ . A* , X* , . J< 2 " 

1 + I - + ~ +••• + *Th 



Inverting the series in we obtain 

*' = 3 ij- - f- (f^f - '-851 (f^)*+ • • • 

__ vP 2 

Note that is an eigen function in irean free path measure and the 

— K ^ I iD y 

function in cm. measure is fc where K s • 



Given an experimental value of K (which is simply the logarithmic 
slope of the flux distribution plotted against 7i ) one may choose a value 
of Z , compute ^ by<K=/^, and using the series in (or any conveniently 

obtained solution) compute « Returning to the integral equation 



V&) ■ 2X / f(l') E '(| 2 '-z|) di ' 



one may insert the computed value of • It is now required to com- 

pute the eigenfunction (ffe} corresponding to the selected value of £ • 

A choice of method is available 0 

LeCaine (6) has treated the problem by a variational method, presenting 
tabulated values for the flux normalized to unit flux at the black plane 
for various values of 2.*/ • A rather rough solution for the parameters 

is available by using these results, and indicates that 2. lies between 
2 o 0 cm”* and 2„4 cm”* . 

It was then decided to attempt solution by a direct numerical method, 
which was carried out to completion. In theory the method goes as follows: 
Given the integral equation as before 

<%) = fbf 

■/d 

The integral is replaced by a Riemann sum equation 

^( z l) = ^( 2 n) (|£.-?J) 

■ h 

which may be put in matrix notation 
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Zo 




E 00 E 01 E„ * - - Eo„ 


z, 




E, 0 E„ E,j - - - - 


* 


- A2 

21 


^20 ^ 2 1 ^22 - - - - 


A 




E Lo ” " ^cn 



where is the absolute value of the exponential integral of argument 

(t- ri) AH .At the argument i = n the exponential integral is singular 
so a mean value was computed for this range. One requires: 




= O'*'**) + E '(f.) 



The integrating scheme is the equivalent of the trapezoidal rule but 
employing the function of mean argument instead of the chordal mean as 
the summation element. 

One proceeds by inserting a tabular trial function, performing the 
matrix operation as indicated and emerging with a new and bettered function 
which is iterated until the variation produced by iteration is insignificant 
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In practice minor departures from the straight forward process were intro- 
duced to reduce the labor. The principal one of these consisted in iterating 
for only about 10 terms of the flux vector and writing the complete vector 
by smoothing the first differences. 

A variant of this scheme was considered, which involved writing the 
equation in a perturbation forms 



^( 2 ) 5 (*) ” 221 5 (*)] dz 

-'-oo 

where is the asymptotic solution, $ is the boundary depression, and £ 
is the perturbation of the operator E t required to extend the integration 
over all 2 . 



= ijrj {[ E '(ii-zi)~ ? (i*'-zi)] 5 (?) +£( ^} d * 

-OO 

or in Riemann sum form 

The equivalent matrices, all infinite order ares 



_ 








- - 






P-2 




E.j-1 ^-21 ^- 2 Z 


v., 


8= 


P-, 


£ = 


- E., 2 E. m E.„ E.,2 - 


p. 

& 




5, 




- 0 0 0 
■ E 7 _, 0 0 0 


- 








- - 0 o 0 


. . 








- 
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E-£ 



000 00000 
oooooooo 
oooooooo 
oooooooo 

= o 0 0 E„ E„ E 1S - - 
o o 0 Ez, E* E i5 - - 

0 0 0 E 5i Ej 2 " - 

0 0 0 - - - - - 

0 0 0 - - - - - 



The first method has the disadvantage that the term products converge 
rather slowly for increasing £ so that each term of the iterated vector 
represents very many term products and the computation is quite laborious. 
The second method avoids this trouble since S vanishes strongly with in- 
creasing 2 , but 5 varies so rapidly near (£= 0) that trapezoidal inte- 

gration is valid only for a very small interval A? . On balance, method 
one appeared most feasible and was chosen. A value of 2.9 cm * was obtained 
for £ with a very fine fit of the solution to the experimental data. 

To investigate the validity of the integration, the process was then 
repeated with an improved kernel. In this case the matrix elements were 
obtained by computing mean values for the exponential integral with a unit 
weight function for each Al n . This process amounts to correcting for 
variation of the kernel over the interval but ignoring the variation of the 
function. A value of I s 2.5 cm -1 resulted. 

A further improvement in the kernel was then tried. One assumes that 
the function may be linearized with small error over the interval Ae 
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The matrix equation is then written: 



where 2 h is the lower boundary of the n interval, Ai is the width of 
the interval, A^n is the difference £ fife ) — * and 

f is the difference M"] . The equation may be separated into 

two terms as follows: 

%) = Tl" £[%")/ E ' (lZ L -?n+?l) d? 

+ il e ^ r?E '( l z i -z^?i) df ] 



The first integral (letting y = |Z£-Z n + f| = |K h + f|) 
/^Kn+i 




The second integral, dropping absolute values: 

Kn+i 

(y-K n )E, (y) dy 

= -^ KB *'( l + K J+C 1 E l(Wi) -K„ 2 E 1(K>i) ]-K n I 1; „ 

I, constitutes the terms of an operator which operates on the function (ftr\ 
and I 2 an operator which operates on the relative increment of the function. 




The process yielded again an excellent fit to the data but for a 
value Z = 2.3 cm“^. Since the computing had become impractically tedious, 
further work on this line was dropped. 

Frankel and Goldberg (3) list a solution for the boundary flux as 
follows: 

%) = _l ) 



where the asymptotic flux goes as the hyperbolic sine vanishing at an 
"extrapolated" boundary thus? 

W»( H ) = slnh X(2+2 0 ) 

The extrapolation distance is given by: 

Z ° * ^ hFF ^Xlt) dS 

The characteristic equation of the system relates & to ^ , so the 

boundary flux may be expressed as a function of it alone: 

- #Fr' 



l <" 


-r) = 


1 

1 +HP& 












1 

1+2^5 




it 


2t 


. it , 


2/" = 


i ' - 


— + 


15 


+ 35 


‘(2r.+lX2n-0 


' H-2^ 


\o 2- i 




. tf 6 , 


_ 


3W& 


A + 


5 


+ — t. 


' ’ * 2n+l 


1 + 2 ipi 
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The procedure is as follows: 

1. The extrapolation distance is computed by Simpson's Rule for 
a few values of K in the neighborhood of 0.15 - 0.20 as in- 
dicated by the numerical solution, and plotted against Z . 

2. Selecting a value of and the corresponding Z Q , the 
experimental boundary flux is normalized to go as 

sinh (Z + Z 0 ) at large argument. This new normali- 
zation is obtained by halving that flux which was normalized 

$ 2 $ ^ 

to go asymptotically as and multiplying by fc 0 

The result is (pp) . 

3. Insert the (fi(p) in the series equation and solve by trial 
for . 

4. Using the computed values of , return to stage (2) and 
iterate the process. Convergence is rapid. 

The described process yielded £ = 2.10 dt 0.05 cm - * . 
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CHAPTER III 



EXPERIMENTAL PROCEDURE 

The medium in which the experiment was conducted was the approxi- 
mately 6' x 6* x 6’ water tank which is mounted on the thermal column of 
the Oak Ridge Graphite Reactor. A current of thermal neutrons passes 
upward through the tank after filtering through the graphite column below. 
The water content is ion-exchange demineralized, and its nuclear properties 
are substantially those of distilled water. 

A disk of aluminum, faced with .020 inch cadmium, and approximately 
10" in diameter, carried on aluminum legs about 10" long, was lowered with 
a rod to a reproducible central location on the tank bottom. A foil holder 
of aluminum and lucite, fastened to the cadmium face was provided with 
adjustable mountings for positioning foils along the central normal to the 
plane. A foil was also cemented to the center of the cadmium plane for 
each exposure. The foil mountings were of lucite of minimum section con- 
sistent with stiffness, and a film of paraffin was employed as an adhesive 
to mount the foils, thus introducing no non-hydrogenous material in the 
immediate vicinity of the foils. The equipment is shown as Fig. (3). 

Indium foils were used, weighing 25 milli-grams nominal, cut from 
50 mg/crn^ stock. These were cleaned of the adhesive wax after each use, 
weighed and allowed to decay for at least 48 hours before re-use. 

Foil counting was done with an end mica-window Geiger-Mueller counter, 
connected through a plate-loaded pre-amplifier stage (of gain ^ 12) to a 
scale of 64. Counter and pre-amplifier were enclosed in a 2" lead shield, 
which reduced background to a negligible fraction of the foil activity. 
Sufficient counts were obtained in each counting run to reduce statistical 
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fluctuation to a negligible level. A number of determinations of counter 
resolving time were made by following the decay of 54 min indium 116 over 
about 5 hours. An average value of 240 yu.sec. was obtained by solving 
the followings 

A, » Ai 

I -N,r I -NjT 

where A = extrapolated activity 

N - mean counting rate, sec - * 

7- z resolving time, sec. 

A correction factor l/l-NT- was plotted against counting rate and 
applied to all counting runs to correct for coincidence losses. 

For exposure, the foil holders (after mounting the foils) were adjusted 
approximately to position with respect to the cadmium plane using a common 
scale. A vernier height gauge was then used to determine the exact foil 
position (estimated probable error ±. .05 mm.). The plane and foil holder 
assembly was then lowered quickly to the tank bottom and exposed for times 
ranging from 5 to 15 minutes depending on foil position to obtain reasonable 
counting rates. No attempt was made to measure either exposure or foil 
’'cooling” times accurately since one foil always and two foils usually were 
exposed at fixed positions to monitor the activity of the variable foil. 

The positions at x = 0 and x= 5.18 cm. are the monitor positions and data 
at these points represents an average of about 15 measurements. In addition 
to the exposures made with axial foils, two exposures were made with foils 
set in a plane parallel to the cadmium to demonstrate that the geometry was 
actually one dimensional. A plot of the diametral variation of flux appears 
as Fig. (6). The maximum deviation from the axial flux amounts to about 
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5% reduction at 8 cm. radius. This is considered satisfactory. To test 
the flux for fast components, a cadmium difference on Indium was taken. 
The ratio of activity in Cd. to activity in the "shadow" of the cadmium 
envelope (a most conservative condition) exceeded 5,000 - therefore foil 
activation due to the 1.44 volt resonance may be safely neglected. 

To analyze the data by the numerical method, it was first plotted 
against the coordinate x in cms. on semi-log paper and the fair curve 
drawn by eye. A clearly exponential branch was obtained beginning at 
x '•'*4.8 cm. The decay constant k of the exponential was solved by 




and fixed at k = 0.375. The scale of the values of (^(x) was then adjusted 
so that (P (asymptotic) = fc °<>3 75x # a trial value of X was then selected, 

and the data plotted, faired and tabulated as a function of 2 in mean path 
measure. The eigenvalue was computed, the tabulated 

on a strip for Az = .2, and numerical iteration of the equation carried 
out until it was obvious that the solution near the plane at E ; 0 would 
not fit the data. The trial was then terminated and another trial made for 
a different value of X . 

As the solutions were improved it became clear that a correction was 
necessary for the remainder contributed by the terms beyond the last tabular 
value of the Eq . The correction was inserted by computing the remainder 
with respect to a proper exponential function and dividing by the value of 
this function at the first untabulated term of the E^ . This number is 
then used to terminate the E^ tabulation. Since the last tabulated value 
of the E^ was for *mAZ = 7.4 one obtains; 
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which may be integrated by writing the E]_ as an integral, then inverting 
the order of integration to yield 



Iteration was carried out by computing only every fifth term of the 
trial function s then interpolating the new function by making the first 
differences regular. The final integration was done in detail in the 
region close to the black boundary. 

Plots of the flux as a function of distance from the black boundary 
appear as Fig. ‘4) and (5). The curves were obtained by using the Le Caine 
results to interpolate the data near the boundary, and the numerically 
iterated function in the region more than 0.5 cm. from the boundary. The 
latter process is considered reasonable since the solutions in this region 
were invariant to changes in the integrating scheme, they are constrained 
by the approach to the asymptotic function, and they fit the data quite well 
Desirable extensions to the analysis might include: 



(1) Computation of the Wiener-Hopf solution for perhaps 15 or 
20 points to firm up confidence in the fit. 

(2) Extensive analysis of the numerical iteration method to estab 
lish the level of computing cfeiail necessary to yield accurate 
solutions using the true values of the parameters. 



Both propositions require computing machine service to reduce program time 
to a practical value. 




which must be computed individually for each value of £ . 
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FIG. I 

VELOCITY DIAGRAM FOR TWO BODY COLLISION 
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FIG. 2 

INTEGRATION FOR PLANE KERNEL - TRANSPORT EQUATION 
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FIG. 3 

FOIL MOUNT AND CADMIUM PLANE ASSEMBLY 
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COMPUTATION OF ACTIVITIES 
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APPENDIX II 



The flux is tabulated as a function of for two values of the 
total cross section. These are solutions corresponding (for *£. s 2,3) 
to the most elaborate numerical treatment, and (for %. = 2.1) to the best 
reconciliation of Wiener -Hopf , variational, and numerical and experi- 
mental solutions. 
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